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We explicitly calculate the moments tn of general Heisenberg Hamiltonians up to eighth order. 
They have the form of finite sums of products of two factors. The first factor is represented by 
a (multi-)graph which has to be evaluated for each particular system under consideration. The 
second factors are well-known universal polynomials in the variable s(s -f 1), where s denotes the 
individual spin quantum number. From these moments we determine the corresponding coefficients 
of the high-temperature expansion of the free energy and the zero field susceptibility by a new 
method. These coefficients can be written in a form which makes explicit their extensive character. 
Our results represent a general tool to calculate eighth-order high-temperature series for arbitrary 
Heisenberg models. The results are applied to concrete systems, namely to magnetic molecules with 
the geometry of the icosidodecahedron, to frustrated square lattices, and to the pyrochlore magnets. 
By comparison with other methods that have been recently applied to these systems, we find that 
the typical susceptibility maximum of the spin-s Heisenberg antiferromagnet is well described by 
the eighth-order high-temperature series. 



I. INTRODUCTION 

The Heisenberg model 

is the basie model to describe physical properties of mag- 
netic insulators. Despite its simplicity the thermodynam- 
ics of the model is generally unknown. For unfrustrated 
quantum spin systems the quantum Monte Carlo (QMC) 
method provides accurate numerical results for the tem- 
perature dependence of the physical quantities. If the ex- 
change couplings are frustrated the so-called "sign prob- 
lem" precludes accurate QMC calculations For one- 
dimensional (ID) frustrated systems the density-matrix 
renormalization group approach^ yields precise results 
in the whole temperature range. For frustrated quan- 
tum spin systems in dimension D > 1 accurate methods 
to calculate thermodynamic properties are notoriously 
rare. Quite reasonable results for arbitrary temperatures 
T can be obtained, e.g.. by a second-order Green func- 
tion technique, see, e. g. , Ref. I^-HI- However, the ap- 
plication of this method needs quite a lot of technical 
experience. Hence, a simple but universal approach is 
desirable. A well established method fulfilling this crite- 
rion is the high-temperature expansion (HTE). Since of- 
ten experimental results, e. g. , for the susceptibility, are 
available in a wide temperature range (including tem- 
peratures exceeding the energy scale set by the major 
exchange constant J, i. e. , for kT ^ |J|), the HTE can 
serve as a method to extract the exchange constants of 
the Heisenberg model from experimental data. 
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For Heisenberg models on the simple two-dimensional 
(2D) and three-dimensional (3D) lattices the HTE is 
available up to high orders, see Refs. and and ref- 
erences therein. However, often one is faced with materi- 
als where two or even more different exchange constants 
are relevant. A typical example are frustrated quasi-lD 
or quasi-2D magnets where except the nearest-neighbor 
(NN) and next-nearest-neighbor (NNN) in-chain or in- 
plane couplings also the interchain or interplane cou- 
plings are important. Typically, for such more complex 
exchange geometries the HTE is known only up to low 
order. In this situation it would be desirable to have 
at one's disposal explicit formulas of higher order HTE 
for general Heisenberg systems and general spin quantum 
number s. It is the aim of the present paper to derive 
such formulas. The key notion is given by the the mo- 
ments Tr " of order n, which can be expressed as sums 
over suitable sets of graphs. From the moments one can 
derive the coefficients of the HTE for, say, susceptibility 
or specific heat in a tedious but straightforward man- 
ner. Unfortunately, the number of involved graphs grows 
super-exponentially with the order n, which delimits the 
maximal order of the HTE for practical purposes. In this 
paper, we have confined ourselves to calculations up to 
eighth order and have to take account of 1139 relevant 
graphs. Nevertheless, this order is sufficient to describe 
typical properties of frustrated spin systems, as we will 
show by means of examples. 

The calculation of the HTE for spin systems has a long 
tradition. Since the 1970s it is known that the moments 
of certain spin lattices with only one exchange constant 
can be written as sums over sets of graphs Qi, with two 
factors. The first factor was called the "lattice constant" 
and counts how often the graph Q^, can be embedded 
into the spin lattice. The second factor is a universal 
polynomial p^{r) in the variable r = s{s + 1). The poly- 
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nomials p^{r) up to eighth order together with the cor- 
responding graphs Gi, are contained in the appendix of 
Ref. @. We have independently calculated these poly- 
nomials by computer-algebraic means and confirmed a 
sample of the data in Ref. 0. The generalization of these 
results from simple spin lattices to arbitrary Heisenberg 
models is achieved by replacing the above-mentioned "lat- 
tice constant" by an "evaluation" of the graph Qi, for 
the spin system under consideration. This evaluation in- 
volves sums of products of coupling constants J^,y and 
yields analytical expressions for the moments of H and 
the coefficients of the HTE of susceptibility and specific 
heat. It seems that such general analytical expressions 
for moments have only be published up to order three, 
see Ref. Those papers that consider higher-order ex- 
pansions, see, e. g. , Refs. [9l-[27j are usually confined to 
special cases, i. e. , special geometries or special values of 
s. We have used some of these papers, namely Refs.l9l.[l0l 
and [2^ to check our general results. 

The paper is organized as follows. In Sec. |ll] we give 
the definitions used and illustrate the underlying mathe- 
matics. In Sec, mil we present general results of the HTE 
coefficients up to fourth order for the moments of the 
Hamiltonian, the free energy, the specific heat, the mag- 
netic moments, and the susceptibility. The very general 
expressions up to eighth order can be found in Supple- 
mentary Material l[2i| of Ref. 0. In Sec. |TV]we apply 
our method to specific Heisenberg models, which are cur- 
rently discussed in the literature, namely the Heisenberg 
antiferromagnet on the Archimedean icosidodecahedron, 
frustrated square-lattice Heisenberg model as well as the 
Heisenberg model on the pyrochlore lattice. For these 
models the HTE for the specific heat and the susceptibil- 
ity up to eighth order for arbitrary spin quantum number 
s are collected in the appendices and Supplementary Ma- 
terial 2113 of Ref. m. 

Although, the information provided in this paper and 
the supplementary materials[2g allows, in principle, to 
calculate the HTE up to eighth order, it might be a te- 
dious task to do so in practice. Hence, we provide a 
simple computer program written in C-| — h that allows to 
calculate within a few seconds the eighth-order HTE co- 
efficients as well as the Fade approximants for the suscep- 
tibility and the specific heat for an arbitrary Heisenberg 
model with up to four different exchange constants (Slj. 



where S denotes 



of H arc defined by fi^ — — ,^ .. , jy — 
the total spin vector and S*^*-*, i = 1, 2, 3, its ith compo- 
nent. As usual, xW) = P TK-(cxp('-ffg)f denotes the 
normalized zero field susceptibility. x{P) = Y^^=i (^nl3^ 
is its HTE in terms of the dimensionless inverse tem- 
perature /3 = jpj^, where J is a typical energy. The 
Hamiltonian H is understood to be dimensionless upon 
division by \J\. The free energy F{/3) is defined by 
-PF{I3) = In (Tr e'^") and its HTE is given by 
—l3F{j3) = S^^o^"!^"- From this one derives the 
normalized specific heat C(/3) = —2/3^^ — and a 

short calculation shows that its HTE C(/3) — Yl^=2 '^nP^ 
is related to that of F(/3) by (i„ = n{n — l)a„ for 
n = 2,3,.... 

HTE are usually written in a compact way by utilizing 
graph-theoretic notations, see, e. g. , Refs. [land 0. Let 
Q he a multigraph consisting of g nodes (vertices) and a 
number of M{i,j) = ^/{jji) bonds (edges) between the 
ith and the jth node. We do not consider "loops," i. e. , 
Af{i,i) = for all i = 1, . . . , 5. The total number of all 
bonds, 7(f/) = J2i<j ^i^^j) ^^^^ ^e called the size of Q. 
Q is not necessarily connected, see the examples below. 
We will identify the set of g nodes with {1,2,..., 17} 
and the set of N spins with {1,2,..., N}. To simplify 
the wording we will omit the prefix "multi-" and simply 
speak of "graphs" in what follows. A selection of graphs 
Q^, V = 1, . . . , needed for purposes of illustration is 
represented in table H] A complete list of all relevant 
graphs up to size 8 can be found in Supplement l.[28| 

For every graph we define its multinomial factor by 



(2) 



Define the symmetry group G{Q) of a graph in the 
obvious way 

G{g) = {7re5g|AA(*,j) =AA(^(*),^(.7)) 

ioT alll <ij < g} . (3) 

Here Sg denotes the group of all permutations tt : 
{l,...,g} — > {l,...,g}. A localization of a graph Q 
is an embedding 



II. DEFINITIONS 



In this paper we consider systems of A'^ spins with 
individual spin quantum number ,s ~ i,l,|,.... The 
Heisenberg Hamiltonian has the form (HI) where the 
JfiL' = Jv^n 1 £ M 7^ < A'^ are suitable coupling 
constants and denotes the spin vector operator of 
the /ith spin. The moments of H will be normalized 
by division by the dimension of the total Hilbert space, 
i. e. 



_ Tr(g" 
- (2s+l)" 



. Analogously, the magnetic moments 



j:{l,...,g}^{l,...,Ar} 



(4) 



up to symmetries of Q. More precisely, two embeddings 
]ii]2 ■ {1, ■ • ■ , ff} — > {I7 • ■ • > N} are called equivalent 
if and only if ji ~ 32 ° for some tt G G{g), and a 
localization of is a corresponding equivalence class 
of embeddings. The number of localizations of Q (for 
given N) will be denoted by L. We will also speak of 
localized graphs Q which will be represented by attaching 
numbers of different spin sites to the nodes of Q, with 
the understanding that two localized graphs which only 
differ by a symmetry permutation of the spin sites are 
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considered as identical, e. g. 



1 



1 



Two localized graphs Qi , Q2 can be soldered in a natu- 
ral way yielding the "soldering product" Qi®Q2, which is 
another localized graph. The nodes of Qi © Q2 are iden- 
tified according to their numbering and the bonds are 
correspondingly added. For example, 



(5) 



Conversely, we will say that the localized graph Qi © Q2 
is decomposed into Qi and Q2- In general, a localized 
graph can be decomposed into different ways. 



From the expansion 



Tri7" = 



Tr 



(6) 



lll<iyi,...,fJ,„<l^n 



it is clear that the expressions for the moments tn involve 
various products of coupling constants J^^. The struc- 
ture of these products can be represented by the graphs 
Q defined above, such that the factors Jf^^ correspond to 
the bonds of Q with multiplicity i. The sum of different 
products in ([6]) of the same structure will be obtained 
by an evaluation of Q, denoted by Q, for the spin sys- 
tem under consideration. Q denotes a real number which 
depends on the coupling constants and only implicitly 
on the number N of spins. This number will be defined 
according to the following statements: 

1. If g > TV we set ^ = 0. 



2. li g < N we select from each equivalence class of 
embeddings a certain representative 

ji:{l,...,g}^{l,...,N}, e^l,...,L (7) 

and define 

L 



e=l l<i<j<g 



(8) 



Obviously, the definition of Q does not depend on 
the choice of representatives ]i since the product 



<i<j<g {'^Je{i),JeU) 



is invariant under permuta- 



tions from the symmetry group tt £ G{Q). 

In order to illustrate this definition we consider an ex- 
ample of = 4 spins and Q = ■ , hence g = 3 < 4 = TV. 
The symmetry group G{Q) consists of all permutations of 
{1,2,3} hence \G{Q)\ = 3! = 6. There are 4! embeddings 
J : {1,2,3} — > {1,2,3,4} and L = || = 4 equivalence 



classes from which we choose the representatives 



Ji = (1 

32 = (1 

J3 = (1 

J4 = (1 



1, 2 
1, 2 
1, 2 
2 2 



2, 3^3), 
2,3^4), 
3,3^4), 

3, 3 ^ 4) . 



(9) 
(10) 

(11) 
(12) 



Hence — Jl2J2zJl3 + J23J3iJ2'i''rJd.iJuJlZ + JliJl2J2i- 

The coefficients c„ of the susceptibility's HTE (and 
similarly the a„ of the free energy HTE) will contain 
products of evaluations Qi, Q^^. These expressions can 
be simplified using rules which transform such prod- 
ucts into linear combinations of other evaluations. To 
give an example, we consider Qi Q2 = ■ . . = 

{^^,<uJ^^y^ {J2k.<\Jk\)- It is obvious that this prod- 
uct can be written as a sum over evaluations of the three 

graphs which can be combined from . . and . . , 

namely . . . . , . . . and . . . In fact. 



(13) 



Similar expressions can be derived for other products 
of evaluations yielding various "product rules" of the form 



(14) 



Here the sum over A runs through all graphs Q\ whose 
localizations are soldering products of localizations of 
and Q^. The integers c^^ count the number of ways to 
decompose a locaHzation of Qx into localizations of 
and Qu- For example, the decomposition ([5]) is unique 
(up to symmetries), hence = 1, c. f. table HI On the 
other hand. 




= 2 



(15) 
(16) 



hence c^q 



2, c. f. table H 

In the case = Qy we have to define c^^ in such a 
way that the binomial factor 2 is included for products 
of different localizations. For example. 



■2 



2 := 



(17) 



From the product rules (|14p one can derive further ones 
for multiple products. 
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TABLE I: A selection of graphs CJ^ 
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14 




16 
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i 1 
1 1 


27 




35 




39 





important rule holds in the case where Q is the disjoint 
union of two simpler graphs, 5 = 5i 1+J ^2 and reads 

v{9)^v(.Qi)v{Q2) .Jf].^ . . (19) 

Note that the polynomials in Ref. @ are defined as our 
•pi, divided by the multinomial factor ([5]), hence these 
factors do not occur in the rule analogous to ([T9l) . Other 
rules, which we need not repeat here, say that the 
vanish a 'priori for certain graphs. 

For the determination of tn it thus suffices to enumer- 
ate the graphs Qv, v G Tn and the corresponding poly- 
nomials Pij. We will give the first four moments for the 
sake of illustration and defer the lengthy expressions for 
t„, n = 5, 6, 7, 8 to the Supplemental Material 1. j28| 

h = (20) 

^2 = = ^^^^ 

h = -^r^^+^r^e^, (22) 

U = 4^'(2-2r + 3r2) g^+^r3(-l + 3r) 
15 9 

^^^^ 



III. RESULTS 
A. Moments 

It turns out that the moments tn can be written in the 
following way: 

tn^Y. '^-P'^i^) ■ (18) 

Here the Qi,, i' G Tn, denote certain graphs of size 
n and the pi, are polynomials of order < n in the 
variable r = s{s + 1). Actually, the pu are of the form 
Pi, = ^i'^^ where g denotes the order of and 

some a^"'* may vanish. The leading coefficients ai"^ 
determine the classical limit r — > 00 of the moments, 
hence they can be calculated by means of integrals over 
unit spheres. 

It is crucial that the polynomials p^, depend neither on 
N nor on the coupling constants J^j^ whereas the terms 
Gd depend only on the coupling constants and only im- 
plicitly on TV via ([5]). The polynomials Pi, up to eighth 
order are well known and have been used for the HTE of 
certain spin lattices. A subset of the p^, is, for example, 
listed in Ref. Q together with certain rules which permit 
the calculation of the remaining polynomials. The most 



B. Free energy 

It is well known that the coefficients of the power series 
for the free energy F{/3) 

00 

- l3Fil3) = In (Tr e"^^) = ^ a„/3" (24) 

n=0 

can be expressed in terms of the moments tn and its 
products. As indicated in Sec. HIl a variety of product 
rules can be used to simplify the resulting expressions. 
This simplification, which is sometimes also referred to 
as the "cumulant expansion", see, e. g. , Ref. 0, has the 
further advantage that it reveals the extensive character 
of the On. By this we mean the following. If the spin 
system under consideration would have a periodic lattice 
structure of, say, K unit cells with periodic boundary 
conditions, it follows immediately that the evaluation of 
a single graph Q linearly scales with K, and hence with 
N, as long as G is connected. For unconnected G the eval- 
uation scales with A'"^ where c is the number of connected 
components of G- Obviously, products of evaluations of 
connected graphs Gi^ Gn would scale with . It turns 
out that the elimination of these and higher products 
in the expression for the a„ by means of the rules 
also eliminates the evaluation terms of unconnected 
graphs. This has to be expected on physical grounds, 
since the total free energy of a spin lattice should be 
an extensive quantity, i. e. , linearly scale with K . But 
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it is an additional consistency test of our results that 
the non-extensive contributions to the a„ actually cancel. 



and give the remaining /i„, n = 4, 5, 6, 7 in Supplemental 
Material I.IH 
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(29) 



The a„, ri = 5,6,7 are given in Supplemental Material 

i.m 



C. Magnetic moments and susceptibility 

To obtain the magnetic moments we will adopt a 
special method which is available if one knows the mo- 
ments tn for all values of the coupling constants J^,^. We 
replace H by the one parameter family of Hamiltonians 
Ha = i/+§(S'^ — A^r) . Equivalently we can substitute 
Jjiv ^ J + o. for all coupling constants. The mag- 
netic moments then result from differentiating Tr(772+^) 
w. r. t. a and finally setting a = 0: 



n + 1 



Tr(iJ^'(S2 -iVr)) (30) 



(n -1-1) (2s 4-1) 



N 



(32) 
(33) 
(34) 



2 n 

ir2(iVr-l)g^+^r3^, 

(8 - (8 + 57V)r + 12r2) ^ + ir3(_l + 2r) Ql 
lr\-\^NT)gl^^r^g-^^'^-v^Q^^. (35) 



The coefficients of the high temperature expansion of 

X = Trr(cxp('-;8ff))^^^ '^^^ expressed through the ^„ 
and the i„ which occur as coefficients of the series in the 
numerator or in the denominator, respectively. The first 
four coefficients are given by: 



n=l 



= /3 - Ml /^^ + 2 (M2 - Mo t2)f 

+ ^(isMo + 3t2Mi ~ M3)/3'* -t- . . . . 
6 



(36) 



(3^„ - Nrtn) 



Inserting the known values for the i„ and the [in yields 
the desired results for the c„. Similarly as in Sec. IIIIBl 
a variety of product rules can be used to simplify the 
resulting expressions revealing the extensive character of 
the c„. 

We will represent the results for the susceptibility's 
HTE up to fourth order in the inverse temperature /?. 
The higher coefficients Cp, n = 5, 6, 7, 8 are given in Sup- 
plemental Material l.[28l| 



We can calculate the left hand side of (pO|) if we in- 
sert the results for the moments and consider "deriva- 
tives" Q' of graphs defined in the following way. Let 
Q'<^^^ denote the graph Q but with one bond removed, 
M{i,j)^ M(i,i)~\. IfA/'(i, j) = Othenwesetg?('^) 0. 
Further let G{Q) and G'(t/(*^^) denote the respective sym- 
metry groups. Then we define 



a' = ^A/-(*,j)e 



\G{Q)\ 



(31) 



ci 

C2 
C3 
C4 



iVr 

2 27^ 

-J Gi, 



9 
1 

18 



02 + ^7-3 "03, 



2 1 

135 ^ T -r ; ^5 ^6 

27 81 



(37) 
(38) 
(39) 

(40) 



One has, so to speak, to break each bond of the graph 
and to sum over all results. Further, one has to introduce 
factors which compensate for the possible change of sym- 
metries. For example, • ' = 6 ■ + . . . . It 

is obvious that the evaluation of Q' just yields ^ ^L-q- 
Then it is a straightforward task to calculate the mag- 
netic moments /iq , ■ • • , Mr by using the above results for 
the tn- We will display the results for /i„, n = 0,1,2,3 



IV. APPLICATION TO FRUSTRATED 
HEISENBERG SYSTEMS 

To improve the HTE approximation G. A. Baker has 
introduced Fade approximants dH (see also Refs. @ 
and 0). These ratios of two polynomials [m, n] = 
Pfn (x) I Rn (x) of degree m and n provide an analytic con- 
tinuation of a function f{x) given by a power series, and. 
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therefore, they yield a better approximation of the func- 
tion f{x). As a rule, approximants with m = n provide 
best results. Since we have a power series up to eighth 
order, we use the corresponding [4, 4] Pade approximant. 



A. The Keplerate magnetic molecules 

In the Keplerate molecules Mo72Fe3o, Mo72Cr3o, 
M072V30, and W72V30 the magnetic ions sit on the ver- 
tices of an almost perfect icosidodccahcdron j33rl36| : see 
Fig. [1] Moreover, the interactions between the magnetic 
ions are well described by the Heisenberg model ([T]) with 
NN interactions. These molecules have attracted much 
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FIG. 2: Susceptibility times temperature in dependence on 
the temperature. The symbols show the experimental data 
for Mo72Cr3o (Ref. [3^ ). the black dashed line represents the 
finite-temperature Lanczos result the blue dashed line 

shows the pure HTE results, and the red solid line the [4, 4] 
Pade approximant of the HTE series. For the exchange pa- 
rameter J and the spectroscopic sphtting factor g we have 
used the same values as in Refs. [3^ and |43. namely J/k = 115 
K and g = 1.95. Note that the Pade approximant and the 
finite-temperature Lanczos data in a wide temperature range 
practically coincide. 



FIG. 1: The Archimedean sohd icosidodecahedron. In the 
magnetic molecules Mo72Fe3o, Mo72Cr3o, M072V30, W72V30 
the magnetic ions occupy the vertices (red bullets). 

attention from the experimental (ssl - issj and theoretical 
side [s^-Hl • One reason is that their frustrated exchange 
geometry has much in common with the kagome lattice, 
see, e. g. , Refs. SI and IH. For the M072V30 and W72V30 
molecules the spin quantum number is s = 1/2, which al- 
lows us to calculate low-energy states exactly by Lanczos 
exact diagonalization [3^. For spin quantum numbers 
s > 1/2 that is already impossible, i. e. , for Mo72Cr3o 
(s = 3/2) and Mo72Fe3o (s = 5/2) the low-energy spec- 
trum can be found only approximately [13, . To eval- 
uate thermodynamic properties already for s = 1/2 one 
has to use approximations [4^ 14^ . Only at high mag- 
netic fields and low temperatures numerical exact results 
were reported [IJ. Very recently a finite-temperature 
Lanczos approximation has been used j4^ to describe the 
magnetic properties of W72V30 at finite temperatures, 
and it has been found that the theoretical results agree 
well with the experimental data over a wide temperature 
range. However, for frustrated quantum spin systems 
with s > 1/2 the calculation of thermodynamic quan- 
tities is even more challenging. Hence our HTE seems 



to be useful in particular for s > 1/2. The HTE series 
for the susceptibility and the specific heat for arbitrary 
spin quantum number are given in Eqs. (jAip and (jA2p 
in Appendix TK] 

We focus on the analysis of the HTE data for the sus- 
ceptibility, since the high temperature magnetic part of 
the specific heat often cannot be accurately separated 
from the phonon part. First we coinpare our s = 1/2 
HTE result for x with experimental [36| and theoretical 
[H] data for W72V30. In Fig. phve show the Tx vs. T 
curve as done in Refs. 13 and l36l. While the raw HTE 
data start to deviate from the experimental ones at about 
T = 115K we find an excellent agreement with the exper- 
imental results and the previous theoretical simulations 
if we use the [4,4] Pade approximant. 

Next we compare in Fig. [3] our results for x for the 
spin quantum numbers s ~ 1/2, 3/2, 5/2 relevant for 
M072V30, W72V30, Mo72Cr3o, and Mo72Fe3o. Again we 
show Tx vs. T, since such a plot is used in many exper- 
imental papers, see, e. g. , Refs. HI, HI, Suggested 
by Eq. (jAip given in Appendix [K\ we use a renormalized 
temperature T / s{s 4- 1). i. c. , we show the dependence 
Tx/s{s -f 1) vs. T/s{s + 1) in Fig. H Obviously the 
curves for different s are very close to each other. From 
Eq. (|Aip it is obvious that with increasing spin quan- 
tum number s in each order of /3 = 1/kT the highest- 
order in r = s(s -|- 1) yields the dominant contribution, 
and therefore the plot Tx/s{s + 1) vs. T/s{s + 1) be- 
comes independent of s for larger values of s. However, 
from both figures [2] and [3] the question arises, whether 
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the Tx/s{s + 1) vs. T/s{s + 1) plot is appropriate to 
detect specific features in in particular at low temper- 
atures. Indeed, the plots in Fig. |4] demonstrate that the 
characteristic low-temperature maximum in x is masked 
in the Tx/s{s + 1) vs. T/s{s + 1) plot. The height and 
the position of the maximum in x clearly depend on s. 
From Fig. 2] it is obvious that its position is shifted to 
lower values of Tx/ s{s + 1) while its height is increasing 
with growing s. 
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FIG. 3: Results of the HTE-Pade approximant for the suscep- 
tibility times temperature in dependence on the temperature 
(arbitrary units) for spin quantum numbers s = 1/2, s — 3/2, 
s = 5/2. 
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FIG. 4; Results of the HTE-Pade approximant for the sus- 
ceptibility in dependence on the temperature for the Kepler- 
ate magnetic molecule for spin quantum numbers s — 1/2, 
s = 3/2, s = 5/2 compared with experimental data for the 
s = 1/2 system W72V30 [1^ (symbols) and finite-temperature 
Lanczos results for s = 1/2 [4J] (blue dashed line). For the 
exchange parameter J and the spectroscopic splitting factor 
g we have used the same values as in Refs.fS^ and 13, namely 
J/k = 115 K and = 1.95. 




FIG. 5: Illustration of the exchange paths for the anisotropic 
frustrated square-lattice model (|41[). 

B. The square-lattice Ji-J[-J2-J2 model 

Next we consider spin systems on infinite lattices. As 
an example we focus on the frequently discussed square- 
lattice Heisenbcrg magnet with NN couplings Ji and frus- 
trating NNN bonds J2, the so-called J1-J2 model. This 
system has attracted a great deal of interest as a model 
system to study quantum phase transitions, see, e. g. , 
the recent publications |47h53| and references therein. The 
HTE for the spin-1/2 J1-J2 model was presented in 
Ref. M- 

The interest in this model is also promoted by a num- 
ber of experimental investigations on magnetic materials 
described reasonably well by the J1-J2 model. However, 
in real materials often one is faced with deviations from 
the ideal J1-J2 model. For instance, in layered vana- 
dium phosphates [13, due to low crystal symmetry 
the bonds along the sides and the diagonals of the square 
can be nonequivalent. Hence, in a realistic spin model for 
these compounds one has to consider two independent 
NN and two independent NNN exchange parameters. 

Therefore we consider here a generalized Ji-J[-J2-J2 
model 

H = Ji ^ Si • Sj + J'l'^ Sj • Sj 

(i,3) {id)' 
+ J2 ^ Sj-Sj+J2 ^ Si-Sj, (41) 

where the NN bonds Ji and J[ as well as the NNN bonds 
J2 and J2 are arranged on the lattice as shown in Fig. [S] 
This model is more appropriate to provide a realistic de- 
scription of frustrated square-lattice materials such as the 
layered vanadium phosphates. 

Based on our general formulas we get the coefficients 
of the high-temperature expansion for the susceptibility 
and the specific heat for the generalized model (|4ip . see 
the Appendix |B] and the Supplemental Material 2 [3^ . 
These formulas contain also interesting limits of coupled 
chain systems [56l - l59j obtained by an appropriate choice 
of the coupling constants. 
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0.12 




0.5 1 1.5 2 2.5 3 3.5 4 



kT/s(s+l) 

FIG. 6: Susceptibility x of the unfrustrated Heisenberg anti- 
ferromagnet [i.e. J[ — Ji — 1, J2 = J2 ~ Q in Eq. (gT])] for 
spin quantum numbers s = 1/2, s = 1, and s = 5/2. The 
data of the pure HTE series and the Pade approximant are 
compared with corresponding QMC (taken from Ref. [s^ and 
ED results. 



0.14 




0.5 1 1.5 2 2.5 3 3.5 4 
kT/J, 



FIG. 7: Comparison of ED and HTE-Pade results for the 
susceptibility of the frustrated model (|4ip for spin quantum 
number s — 1/2. For the exchange parameters we follow 
Ref. [5^ and choose Ji = J{ = — 1/2 and 

{J2 + J2)/{Ji + J\) = —2. The thermodynamic energy scale 
is defined as Jc = ^( + + J| + J2^)/2. 



First, we compare in Fig. [5] the HTE data for the sus- 
ceptibility with accurate QMC data for the pure square 
lattice Heisenbergantiferromagnet for s = 1/2 and s = 1, 
see, e. g. , Refs. |63, HH andlssl as well as with numerical 
exact data for finite lattices obtained by full exact diago- 
nalization (ED). Again we use the renormalized tempera- 
ture T/s{s + 1) for the plot; see the discussion in the pre- 
vious section. The comparison with precise QMC data 
allows to estimate that temperature Ta down to which 
the HTE approximation for x is accurate. We find that 
the pure HTE in eighth order practically coincides with 
QMC data until Tas/{Jis{s + 1)) = 1.2, 1.4 and 2.0 for 
s = 1/2, s = 1, and s = 5/2, respectively. Using the [4,4] 
Pade approximant we find Ta^2/ {Jis{s + 1)) = 0.85, 0.75, 
and 0.95 for s ~ 1/2 s = 1, and s = 5/2, respectively, 
and it is evident from Fig. [S] that the maximum in x is 
described accurately. Even significantly below Ta.2 the 
Pade approximant describes the QMC data reasonably 
well. Moreover, by comparison with ED data we can fig- 
ure out how good typical ED results can describe realistic 
large systems in two dimensions. Often, the ED is used 
as the only method to discuss the thermodynamics of 
strongly frustrated 2D quantum spin systems, see, e. g. , 
Refs. I54 and[63-[63. The results shown in Fig. [5] indicate 
that for 2D systems already at moderate temperatures 
and even for s = 1/2 (where largest systems are accessi- 
ble by ED) significant finite-size effects appear, and that 
our HTE results for N ^ 00 are better then typical ED 
results. A similar finding was reported in Ref. |5| where 
ED results for x are compared with data of a Green's 
function approach for a spin-1/2 frustrated square-lattice 
fcrroniagnet. 



We consider now the generahzed Ji- J{- J2- J^_model 
(|4T|) relevant for layered vanadium phosphates [53, [El] . 
First, we mention that, for the symmetric model (i.e. 
Ji = J[, J2 = J2), we give the general formulas for the 
HTE coefficients for arbitrary s up to eighth order in Ap- 
pendix |B] For the asymmetric model for arbitrary s the 
formulas become very lengthy for higher orders. There- 
fore, in Appendix [B] we present the formulas for arbitrary 
s only up to fifth order, and give the remaining sixth to 
eighth orders in Supplemental Material 2.^^ To illus- 
trate our HTE results we follow the lines of Ref. [13 and 
discuss the influence of exchange asymmetry Ji =^ J{, 
J2 7^ J2 oil the temperature dependence of the suscepti- 
bility, in particular, on the position and the hei ght of the 
maximum in x- This issue was discussed Ref. |5J based 
on ED data for = 16 = 4 x 4 (see Fig. 9 therein). We 
have repeated these ED calculation and compare the ED 
results with the HTE data for = 16 and A^ 00 in 
Fig. Hand in Table [Hi 

Obviously for A^ = 16 the ED and the corresponding 
HTE-Pade data for the maximum in x agree well. But 
it is also obvious, that the finite-size data for the max- 
imum do not agree well with data for N ^ 00. The 
shift of the maximum by varying the asymmetry (i. e. , 
the difference in Ji and J[ or/and in J2 and Jj) dis- 
cussed Ref. 113 is not observed (or is at least much less 
pronounced) in the HTE results for A^ — > 00, cf. Fig. [7] 
and Table [III Hence we argue again that the conclu- 
sions based on finite-temperature ED data for 2D systems 
might be not reliable for large systems. 
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TABLE II: Position Tm = Tmax/Jc and height Xm = 
XmaxJc/{NAg'^ fJ-'s) of the susceptibility maximum of the 
model (|4H) . In the table we compare the numerical ex- 
act ED data for A*' = 16 (superscript ED), the corre- 
sponding HTE-Pade data for TV = 16 (superscript HTEia) 
and the HTE-Pade data for N = oo (superscript HTEao)- 
The thermodynamic energy scale is defined as Jc = 

^(j2+J^2+j|+j,2)/2. 



Jl 




J2 


J'2 


rpED 
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rpHTEie 






HTE^ 
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2 


2 





1.02 


0.1214 


1.02 


0.1214 


0.83 


0.1290 
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1 
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8 
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0.99 


0.1199 


1.01 


0.1197 


0.81 


0.1284 


1 
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1 
2 


4 
3 


2 
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1.00 


0.1181 


1.04 


0.1174 


0.80 


0.1287 


1 
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1 
2 


8 
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6 
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1.02 


0.1170 


1.06 


0.1159 


0.80 


0.1291 


1 
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1 
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1 


1 


1.02 


0.1167 


1.07 


0.1155 


0.80 


0.1291 


4 
7 


3 
7 


1 


1 


1.03 


0.1165 


1.08 


0.1152 


0.80 


0.1288 


2 
3 


1 
3 


1 


1 


1.05 


0.1154 


1.11 


0.1141 


0.83 


0.1276 


4 
5 


1 
5 


1 


1 


1.10 


0.1136 


1.17 


0.1118 


0.89 


0.1254 


-1 





1 


1 


1.15 


0.1115 


1.22 


0.1097 


0.96 


0.1225 



C. The Heisenberg model on the pyrochlore lattice 

As the last example we consider a 3D frustrated spin 
system, namely the Heisenberg model on the pyroehlore 
lattice. In three dimensions the ED is not applica- 
ble to calculate reasonably well thermodynamic prop- 
erties. Moreover, typically there is finite-temperature 
phase transition which needs special analysis of the HTE 
series. The pyrochlore lattice is highly frustrated and 
it has attracted much attention over the last years, see, 
e. g. , Refs. IgsI - IbtI and references therein. To the best 
of our knowledge so far no higher-order HTE has been 
presented. For the classical limit the thermodynamics 
was investigated systematically mainly by classical Monte 
Carlo (MC) simulations; see, e. g. , Refs. [H HI and (z^. 
Due to strong frustration there is no phase transition 
to an ordered low-temperature phase for the pyrochlore 
Heisenberg antiferromagnet. For the quantum model no 
precise data are available at lower temperatures. 

The HTE series for the susceptibility and the specific 
heat for arbitrary spin quantum number s are given in 
Eqs. dCl]) and in Appendix O The plots of the 

Fade approximants for the Heisenberg antiferromagnet 
are shown in Fig. |8] for various values of s. For the 
classical model (s oo) we compare our HTE data 
with MC data calculated in Ref. [6^ see also Ref. [tO- 
Surprisingly, there is an excellent agreement with the 
MC data down to temperatures which are considerably 
below \J\/k. In particular, the fact that there is no 
maximum in the x(^) curve is observed both in MC and 
HTE results. Lowering the quantum number s, i. e. , 
increasing the quantum fiuctuations a low-temperature 
maximum in x{T) emerges. The height Xm of the 
maximum decreases, whereas the position Tm/s{s + 1) 
increases with decreasing of s. 
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FIG. 8: HTE Pade approximant for the susceptibility of py- 
rochlore Heisenberg antiferromagnet for s = 1/2, 1, 3/2, 5/2 
and s ^ oo. The MC data for s ^ oo are taken from Refs. l69l 
andH 



V. CONCLUSIONS 

In this paper we provide general expressions for the 
high-temperature expansion series up to eighth order of 
free energy, the specific heat, and the susceptibility for 
Heisenberg models with arbitrary exchange patterns J^i, 
and spin quantum number s. These formulas can be used 
as a tool to investigate thermodynamic properties of gen- 
eral Heisenberg systems and thus for the interpretation 
of experimental data, especially if other precise meth- 
ods, such as the quantum Monte Carlo method or the 
finite-temperature density matrix renormalization group 
approach, arc not applicable. By comparison with pre- 
cise quantum Monte Carlo results for the susceptibility x 
of the unfrustrated 2D Heisenberg antiferromagnet with 
NN exchange J with s = 1/2, s = 1, s = 5/2 we 
find that the HTE results yield the correct susceptibility 
at high temperatures down up to T/s(s + 1) « \J\/k. 
Using Fade approximants, the accuracy can be extended 
to lower temperatures. In particular, the typical maxi- 
mum in X for the Heisenberg antiferromagnet can be well 
described using the HTE of eighth order. 

We apply our method to frustrated systems, namely to 
frustrated Keplerate magnetic molecules, to a frustrated 
square-lattice Heisenberg magnet, and to a pyrochlore 
Heisenberg magnet. By comparison with finite-size 
data for the unfrustrated as well as the frustrated 
square-lattice Heisenberg model obtained by full exact 
diagonalization we find that the size of 2D systems 
accessible by full exact diagonalization seems to be 
too small to get precise data for the susceptibility 
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maximum. The comparison with Monte-Carlo data 
for the classical pyrochlore Heisenberg antiferromagnet 
yields an excellent agreement down to low temperatures. 
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Appendix B: The high-temperature expansion for 
the susceptibility and the specific heat for the 
square-lattice J\-,J[-J2-J2 model 

Here we list the general formulas for the susceptibility 
and the specific heat for the Ji-J[-J2-J2 model defined 
in Eq. (j4ip . Since the corresponding formulas become 
very lengthy in higher orders of the HTE, we restrict 
ourselves here to (i) general formulas for the symmetric 



Appendix A: The high-temperature expansion for 

the susceptibility and the specific heat for the 
Heisenberg model on the icosidodecahedron with 
nearest-neighbor exchange interaction 



The general formulas for the susceptibility and the spe- 
cific heat for the Heisenberg model on the icosidodeca- 
hedron with the NN exchange constant J up to eighth 
order read 



I 

model ( Ji = J( and J2 = J2) foi' arbitrary spin quantum 
number s up to eighth order and (ii) general formulas for 
the asymmetric model (Ji ^ J[ and J2 ^ J2) for arbi- 
trary s up to fifth order, only (for the remaining sixth 
to eighth order coefficients, see the Supplemental Mate- 
rial 2 [131 )• Note that for the symmetric models with 
s = 1/2 the HTE coefficients up to 10th order are given 
in Ref. ^ 

First we give the formulas for the symmetric model: 



N 



(Al) 



Cl 

C5 

C7 

C8 



3 r ; C2 = --r 



C3 



9 



^(4r-l) ; C4 = - 



4860 



(-45 + 702 r - 1892 + 1328 r^) 



405 



-(3-28r-f 37r2) 



510 300 



2 041 200 

„2 



11 481 750 



(1728 - 35 946 r + 164 289 - 209 296 + 99 776 r^) 

(-2898 -I- 72 972 r - 467 127 + 967 124 - 765 536 + 259 008 r^) 

(7695 - 223 128 r + 1 769 382 - 5 284 101 + 6 231 056 - 3 632 860 + 745 760 r^) 



CX(3)=NkY,dniJPr 
2 



dr 
dn 



71 = 2 



"2=3'- ; ds ^ —{3 - Ar) ; ^4 = 

= — -(9 - 126 r + 116 + 32 r^) 
162 



(A2) 



45 



-3 + 23 r + 3 r^) 



22 680 



97 200 



(576 - 11 142 r + 34 323 + 5088 + 3952 r^) 
(-1242 + 29 556 r - 150 039 + 100 736 + 32 624 + 25 472 



1 093 500 



(-7695 + 213 084 r - 1 435 806 + 2 566 548 -I- 214 682 + 473 600 + 82 120 r^). 
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x(/3) =iV^c„/3" ; r = 5(s + l) 
4 



Cl 
C4 

4 

C6 



n=l 
C2 



9 



r'(Ji + J2) 



C3 



1 

— r 
27 



405 
1 

4860' 



C7 



3061800 



(-4347 + 99738r - 623943r2 + 1392666r3 - 1440944r^ + 673152r^) 

+28 Jf J2r(2133 - 44406r + 223896r2 - 351328r3 + 209664r^) 
+42 Jj* J|r(1620 - 35649r + 248050r2 - 618992r3 + 485984r'*) 
+28 JiV|r(2061 - 46134r + 389520r2 - 1234240r3 + 1254912r^) 
+2 JiV|r(36936 - 789687r + 4517328r2 - 12128272r3 + 15232512r^) 
+64Ji J|r(432 - 11097r + 75108r2 - 181180r3 + leGGSOr^) 

+j|(-4347 + 99738r - 623943r2 + 1392666r3 - 1440944r^ + 673152r5) 
1 



C8 



91854000 



3 (20520 - 5361 12r + 4174761?'^ - 12734370r3 + 18166056r^ 

-13785984r'5 + 5028608r'^) + 8 jfj2r(- 108459 + 2561472r - 17865060r2 + 42056212r3 
-43723408r4 + 19466016r^) + 8 jf j|r(-120123 + 3081636r - 26591049r2 
+104756322r3 - 150488432r* + 82730944r5) + 4 jf j|r(-205416 + 6235740r 
-66481989r2 + 275607666r3 - 543091408r4 + 372472256r'5) + 2 JiV|r(-364932 
+10618461r - 111258000r2 + 530418656r3 - 1116577664r* + 926078848r^) 
+2OJ1 J2V(-51678 + 1313703r - 10031442r2 + 34441608r^ - 66122048r'' 
+61121472r''^) + 80 Ji J2^r(-4347 + 127044r - 1130751?'^ + 3894570r^ - 6116240r^ 
+4059904r^) + 3jJ(20520 - 536112r + 4174761r2 - 12734370r3 + 18166056r'' 

-13785984r^ + 5028608r6) 



and 



C(/3) =iVfc^d„/3" 

2 

5 
2 



d2 
(is 



n=2 
t2 I t2\„2 



1 



3(^1+ ^2 K ; ds = -r\Jt-8JtJ2r + J^: 



45 
1 

I62' 



Jf{3 - 18r + 7r^) - 20jfj2r + 50JiV|r(-l + 2r) + j|(3 - 18r + 7r^) 
Jf (9 - 108r + 28r^) - Sjfj2r{9 - 104r + 16r2) + 40Jf j|r(-3 + 5r) 

4JiV|r(39 - 274r + 256r2) + j|(9 - 108r + 28r^) 
1 



22680 



J{'(-576 + 9756r - 23739r^ + 12160r'' + 640r^) + 560Jf J2r(9 - 135r + 16r^) 



-28Jfj^r{~297 + 4023r - 12904r2 + 824r3) + 280JiV|r(24 - 153r + 128r2) 



(Bl) 



3Ji^(-l + 4r) + 32 Ji Jar + 3J|(-1 + 4r) 

rr- 3Ji^(2 - 17r + 28r^) + lOJ^ J2r(-9 + 34r) + 20Ji j|r(-3 + 20r) + 3J|(2 - 17r + 287'^) 

J]*(-45 + 648r - ISOSr^ + 1712r3) + 120JiV2r(5 - 54r + 80r2) 

12 JiV|r(69 - 574r + 1536r2) + 192Ji j|r(2 - 27r + 68r^) + J|(-45 + 648r - 1808r2 + 1712r3) 

= --TT^^r^ [2 Jf (216 - 4131r + 18339r2 - 28710r3 + 18100r4) + 21 jfj2r(-279 + 4801r - 14048r2 + 12368r3) 
-14JiV|r(-477 + 7158r - 40044r2 + 49864r3) + 35 jf j|r(-216 + 3261r - 13504r2 + 24240r3) 
-350JiJ2ri-9 + 186r - 832r2 + 1168r3) + 2 j|(216 - 4131r + 18339r2 - 28710r^ + I8IOO/) 



(B2) 



12 



-70JiV^r(-144 + 2139r - 3872r2 + 2400r^) + j|(-576 + 9756r - 23739r2 + 12160r3 + 640r 



1 



^ 97200 



Jlil2A2 - 25920r + 116211r2 - 58432r3 - 1920r^) + 8Jf J2r(-1503 + 28233r 
-99398r2 + 31296r3 + 36mr^) + 28Jf j|r(-639 + 11076r - 31508r2 + 1648r3) 

+56Ji J2V(-279 + 5943r - 42779r2 + 59 1 007'^ + 64r^) + 42jj^4r(-306 + 4361r - 7144r2 + 4000r3) 
-4Ji2 J|r(5472 - llOOOZr + 401658r2 - 434720r3 + 186560r^) 

+ Jj(1242 - 25920r + 116211r2 - 58432r3 - 1920r^) 

1 



. J?(7695 - 185976r + 1160352r2 - 1811898r3 + 889724r'' + 20256r^ 
1093500 L 

-27512r^) + 8Ji^ J2r(-10044 + 212661r - 1119348r2 + 371668r3 + 27120r^) 
-2jf J2r(57672 - 1256355r + 7596708r2 - 19368452r'' + 3177376r'* + SWSOSr^) 
+20Jf j|r(-4536 + 111771r - 777870r2 + 988176r3 + 896r'*) - 2 Jj* j|r(44874 
-1369062r + 13508715r2 - 39201672r3 + 31536928r'* + 2049424r^) - 20 j|r(3618 
-69579r + 233613r2 - 237744r3 + 93280r'') + 10 JiV|r(-13662 + 325665r - 1954599r2 
+3232376r=^ - 2405088r4 + 818624r5) + J2^(7695 - 185976r + 1160352r2 - 1811898r=^ 

+889724r4 + 20256r5 - 27512r6) 



Next we give the formulas for the asymmetric model (up to fifth order). For the susceptibility x we find 



(B3) 



^^ = 3 ' 

1 2 

^^=54" 



+AJ[{Ji + 4 
1 



C2 = --iJ'i + Ji + J2 + J2y; 

y 

J! + 4^ + J|) + 4(jf + J2 + + 4j/ + J| + 4,/i(J^ + J2) 



C4 



405 

T/2/ 



+10J['{Ji 
+10 J[{J^ 
+20Jf(Ji 
+5J((4J2 
1 . 



J2))r] ; 

- 6( + + J'i" + 4) + 3(7Jf + 7Jl + 10J2(J^ + J2) 

- + J2) + 10Ji(J^2 _^ j2) ^ _^ J2)(7J^2 + 3J^ J2 + 7J|) 
42 + j| + + j^)))r - 4(jf + Jl + 20J2(J^ + J2) 

- J'2 + J2) + 20Ji(J^2 _^ 3 J/ _^ j2) ^ (J/ j^2 ^ ^gj/ _^ j2) 



C5 



9720 
+12J{r(16Jf 
+87J?(J^ + J2) + 12(J^ 



9Ji(4 + . 
- 45 (Jf 4 
9J2 



4(JJ 
+ J|) 

/2 



3J^ J2 + J|)))r2 
192jf (Jl + + J2)r(l - 6r + 4r2) 



9J1J2 



16J2-'' + 9 J^/2 



9JiJ| + 16 J| - 2(48 



^2)(4J^' 



4J2 + 4J2) + Ji(87jf + 80J^J2 + 87 J|))r 
+16(Ji + J2 + ^2)(4Ji^ + 4J^2 _^ 26J^ J2 + 4J| + 17Ji( + J2)y) ~ J(^(45 + 4r(-69 
+4r(ll + r))) + 12jf r(30(J^2 _^ j2) _ 8o(J^ + J2fr + 80(2J^ + J2)(J^ + 2J2)r^ 
+10Ji2(3 + 8r(-l + 2r)) + 3Ji(J2 + J2)(3 + 2r(-29 + 56r))) + 4r(48J2V2(l - 6r 
+4r2) + 48J2 J|(l - 6r + Ar^) + 48Ji'( J2 + J2)(l - 6r + 4?'^) + 30Ji2(3( jf + J^) 
-8(4 + J2fr + 8(2J^ + J2)( + 2 J2)r2) + j4(69 - 4r(ll + r)) + J2'^(69 
-4r(ll + r)) + J|(69 - 4r(ll + r)) + 30Jf j|(3 + 8r(~l + 2r)) + 24Ji(J^ 

+J2)(2 Jf (1 - 6r + 4r2) + 27^(1 - 6r + 4r2) + J2(-2 + r(-3 + 52r)))) . 
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For the specific heat C we have 
C(/3) 



Nk V 



n=2 



J2 + + j2) 
Jf + J'l + J^^ + 4 - ^J[.h{j'2 + J2) 

3(J{^ + Ji^ + + J|) - 2(4j;* + 2(2 Jf + 5 JiVs^ + 2J^* + b{Jl + Js^) j| + 2 j|) 



+hJ[Ji{J'^ + J| + Ji(J^ + J2)) + 5jf (2J2 + + J2) + 2{j'i + Jl)))r + (-3J{* 

-372"* + 20J2V2 - 3 j| + 20J1 ( J2^ + 3J2 J2 + J2 ) + 20jf ( Ji + j'i + 3J2J2 + Ji)y 

1 



3Jf 



324 

-2J|)2^ 

-3Ji(J2 



9(jf + + J^^ + 4) - 12(4jf + 4jf + 5Jf(J^2 _^ j2) ^ _^ J2)(2J^2 _ J/ 
5J?(jf + J|) + 5J[^{jf + J'i + J|) + 3J( + ^1 + (J2 + -^2)) + J'li^Jl 



J2 



'2 '-'21 
t2\ 



4(-3 Jf - 3 Jf 



3 J'^ 



104V? 



34 



10^(4^ + 34 J2 + J|) + 5Jf (4 + J2)(242 + 4 J2 + 2J|) + 8J( Ji(4 + J2)(13(Jf + 4^) 
~2J'J: 



2^2 



13 J|) 



2jf (542 + 52Ji(4 + J2) + 5(42 _^ 3j/ + + 5Jf (2J3 + (4 



J2)(242 ^ ^ 2j|)))r2 - 64J{ Ji(4 + J2)(jf + j2 + ^ ^ j2)^3 



r 



(B4) 



Appendix C: The high-temperature expansion for 
the susceptibility and the specific heat for the 
Heisenberg model on the pyrochlore lattice 



tice with NN exchange constant J up to eighth order read 
for the susccptibihty as follows: 



The general formulas for the susceptibility and the spe- 
cific heat for the Heisenberg model on the pyrochlore lat- 



J 



Cl 

C4 
C5 

C6 
C7 
C8 



C2 



2r2 
~3~ 



C3 



18 



135 



r2(6 - 91r + 224?'^ 



(Cl) 



r2(-3 + 20r) 



1080 



14175 
1 

2041200' 
1 



r^{-l5 + 376r - ISlGr^ + 2544r3) 

'(72 - 2406r + ISDOOr^ - 47188r3 + 46848r^) 



61236000 
-601876480r^ + 399408640rS) 

and for the specific heat 



^(-4347 + 176346r - 1901709r2 + 7300134r3 - 11982944r'' + 9482624r^) 
(61560 - 2887056r + 38320749r2 - 202461642r3 + 477409712r'' 



(C2) 



4 = Ir^iS - 8r) ; = ^r^3 - 38r + 7r^ 
o 15 



14 



= -^r^(3 - 68r + U8r^ + 32r^) 
So 

de = ---47rT7-^(-1728 + 53964r - SOlGTlr^ + 102672r3 + 56128r'^) 
45360 

d7 = -7ri7;7:7-2(-1242 + 47808r - 418437r2 + 728520r3 + 178240r* + 13312r5) 
64800 

ds = Q r^(7695 - 345816r + 3954204r2 - 13638312r3 + 5728812r'^ + 4024640r'5 
+ 1856680r6). 
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